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Abstract 

We use a coarse-graining approach to prove that inter-scale transfer of kinetic energy in compress- 
ible turbulence is dominated by local interactions. Locality here means that interactions between 
disparate scales decay at least as fast as a power-law function of the scale-disparity ratio. In partic- 
ular, our results preclude transfer of kinetic energy from large-scales directly to dissipation scales, 
such as into shocks, in the limit of high Reynolds number turbulence as is commonly believed. The 
results hold in broad generality, at any Mach number, for any equation of state, and without the 
requirement of homogeneity or isotropy. The assumptions we make in our proofs on the scaling of 
velocity, pressure, and density structure functions are weak and enjoy compelling empirical support. 
Under a stronger assumption on pressure dilatation co-spectrum, we show that mean kinetic and 
internal energy budgets statistically decouple beyond a transitional "conversion" range. Our analysis 
demonstrates the existence of an ensuing inertial scale-range over which mean SGS kinetic energy 
flux becomes constant, independent of scale. Over this inertial range, mean kinetic energy cascades 
locally and in a conservative fashion, despite not being an invariant. 
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1 Introduction 



This paper is the second in a series which investigates the physical nature of compressible turbulence. 
In the first paper [T], we laid a framework to study the coupling of scales in such flows and to analyze 
transfer of kinetic energy between different scales. The aim here is to explore whether this transfer of 
energy takes place through a scale-local cascade process, similar to incompressible turbulence. This is of 
central importance in the subject because scale locality is necessary to warrant the concept of an inertial 
range and to justify the existence of universal statistics of turbulent fluctuations. 

The traditional Richardson-Kolmogorov-Onsager picture of incompressible turbulence makes the fun- 
damental assumption of a scale-local cascade process in which modes all of a comparable scale I 
(differing at most by some fixed ratio, typically of order 2) participate predominantly in the transfer of 
energy across scale L This also implies that energy transfer is primarily between modes at comparable 
scales, with a ratio of order 2. If, furthermore, the cascade steps are chaotic processes then it is expected 
that any "memory" of large-scale particulars of the system, such as geometry and large-scale statistics, 
or the specifics of microscopic dissipation, will be "forgotten." This gives rise to an inertial scale-range 
over which turbulent fiuctuations have universal statistics and the flow evolves under its own internal 
dynamics without direct communication with the largest or smallest scales in the system. 

Therefore, scale locality is crucial to justify the existence of an inertial range and its universal 
statistics, and is necessary for the physical foundation of large-eddy simulation (LES) modelling of 
turbulence. It motivates the belief that models of subscale terms in the equations for large-scales can be 
of general utility, independent of the particulars of turbulent flows under study. 

[29j l30l 131] was the first to demonstrate locality in incompressible Navier-Stokes turbulence using 
detailed closure calculations. He showed that interactions between widely separated scales, ^ ^2, 
decay as a power-law of their ratio, {lijl-i)"' ^ where a > 0. Later on, [1^ was able to prove local- 
ity rigorously from the equations of motion and under very mild assumptions, without any closure or 
statistical averaging. More recently, [2] proved locality using Fourier analysis and in [3 , they showed 
that it also holds in incompressible magnetohydrodynamic turbulence. Furthermore, there has been sev- 
eral recent studies by []Domaradzki07a,Domaradzki09,AluieEyink09,AluieEyinklO,DomaradzkilO which 
support the aforementioned theoretical results from direct numerical simulation (DNS) data. 

No similar results, either theoretical or empirical, exist for compressible turbulence. The idea of 
a cascade itself is without physical basis since kinetic energy is not a global invariant of the inviscid 
dynamics. The notion of an inertial cascade-range in such fiows is, therefore, still tenuous and unsub- 
stantiated. In this paper, we shall prove rigorously under modest assumptions that inter-scale transfer 
of kinetic energy is indeed local in scale. Under a stronger assumption, we will further show that kinetic 
energy cascades conservatively despite not being an invariant. 
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The outline of this paper is as fohows. In §[5]we present prehminary definitions and discussion. In the 
course of proving locality of the cascade, we shall first prove that the subgrid scale fiux (defined below) 
is dominated by scale-local interactions. This is done by expressing the flux in terms of increments 
in §[21 then proving that increments themselves are scale-local in §31 then, flnally, showing how this 
leads to locality of the flux in §[5l In §[6l we discuss the implications on locality of the cascade itself. 
We show that beyond a transitional "conversion" range, an inertial range emerges over which kinetic 
energy cascades conservatively in a scale- local fashion. In §[7jT]and §[721 we argue for the validity of our 
assumptions based on empirical evidence and physical arguments. In ii l7.4l and ^ I7.5l we discuss common 
misconceptions regarding scale locality in the presence of intermittency and shocks. We summarize our 
main results in §[51 and defer mathematical details to Appendices O and |B1 

2 Preliminaries 

2.1 Governing dynamics 

We prove locality of kinetic energy transfer by a direct analysis of the compressible Navier Stokes equa- 
tions, without use of any closure approximation. The equations are those of continuity momentum 
([H, and either internal energy ([3]) or total energy ([4]): 

dtp +d,{pu,) = Q (1) 

dt{pui) +dj{pu.iUj) = -diP + pdj{dju.i + ^dmUmSij) + pFi (2) 

dt{pe) +dj {peuj - p,{umd„iUj ~ UjdmUm)} = ~PdjUj + p\djUi\'^ + ^IdjUj]"^ - djqj (3) 

2 

dt(pE) +dj{pEuj) = ~dj{Pu,j) + pdj{ui[{djUi + diUj) - -dmU,nSij]} - djq-j + puiFi (4) 



Here, u is velocity, p is density, e is internal energy per unit mass, E = |up/2 + e is total energy per 
unit mass, P is pressure, p is dynamic viscosity, F is an external acceleration held stirring the fluid, 
q = —K^T is the heat flux with a conduction coefficient k and temperature T. For convenience, we 
have assumed a zero bulk viscosity even though all our analysis applies to the more general case. We 
have also assumed that p = vp la independent of x. 

2.2 Coarse-grained equations 

Following 25^ and I l9j . we presented in a previous paper, [1 , a scale-decomposition based on coarse- 
graining which satisfies an inviscid criterion, i.e. it guarantees that viscous momentum diffusion and 
kinetic energy dissipation are negligible at large-scales. The decomposition yields a scale-range L 3> ^ 3> 
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£^ over which kinetic energy is immune from viscous dissipation and external injection by stirring. Here, 
L denotes "integral scale" and £^ denotes dissipation scale. 

Using the coarse-graining approach, we can resolve dynamics both in scale and in space. We define 
a coarse-grained or (low-pass) filtered field in d-dimensions as 

a,(x) = J d''rGiir)a{^ + r), (5) 

where G(r) is a smooth convolution kernel which decays sufficiently rapidly for large r, and is normalized, 
J d'^r G{r) — 1. Its dilation Gi{r) = £^'^G{r/£) has its main support in a ball of radius £. We also define 
a complimentary high-pass filter by 

a^(x) = a(x) -a£(x). (6) 

In the rest of our paper, we shall take the liberty of dropping subscript £ whenever there is no risk of 
ambiguity. 

In [1], we proved that viscous momentum diffusion and kinetic energy dissipation are negligible at 
large-scales when a Favre (or density- weighted) decomposition is employed. A Favre filtered field is 

weighted by density as 

/.(X) ^ (7) 

Pern 

The resultant large-scale dynamics for continuity and momentum are, respectively, 

^tp + ^^{pu^)^0. (8) 



dipui + dj{pui Uj) = -dj ijD f{ui,Ujj) - diP 

2 ~ 

-I- p.djilidjU^^' diUj) - -drnUjnSij]} ^Wi^ (9) 

where 

'pf{ui, Uj) = p{UiU-j ~ Ui Ui) (10) 

is the turbulent stress from the eliminated scales < £. It is also straightforward to derive a kinetic energy 
budget for the large-scale, which reads 

9tP^-^ + V-Jf = -n, - A, + PeV-ue - De + e"^ (11) 

where Jf(x) is space transport of large-scale kinetic energy, lle{x.) + is the subgrid scale (SGS) 

kinetic energy fiux to scales < £, — PV-u is large-scale pressure dilatation, Di{x.) is viscous dissipation 
acting on scales > £, and e'"'' (x) is the energy injected due to external stirring. These terms are defined 
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as 



n,(x) 

Af(x) 



p djUl f{ui,Uj) 



-djP t{p,Uj) 



(12) 
(13) 

(14) 

(15) 
(16) 



where 

Mf,g) = UWi-7i9i (17) 

in expression (1131) is the 2"'^-order generalized central moment of fields f{x.),g{x) (see |25jV Equations 
(|8))- (ITT]) describe the dynamics at scales > £, for arbitrary £, at every point x and at every instant in 
time. They hold for each realization of the flow without any statistical averaging. 

The SGS flux is comprised of deformation work, Hi, and baropycnal work, , which we discussed 
in some detail in [T] . These represent the only two processes capable of direct transfer of kinetic energy 
across scales. Pressure dilatation, — P^V-Uf, does not contain any modes at scales < i (or a moderate 
multiple thereof), at least for a filter kernel G'(k) compactly supported in Fourier space. Therefore, 
pressure dilatation cannot participate in the inter-scale transfer of kinetic energy and only contributes 
to conversion of large-scale kinetic energy into internal energy. This observation is one of the key 
ingredients to proving scale locality. 

3 SGS flux in terms of increments 

As was realized in the pioneering work of [19], there are two facts crucial for scale locality of the SGS 
flux across i. First is that SGS flux can be written in terms of increments, 

5a(x;r) = a(x + r) - a(x), (18) 

for separation distances |r| < i (or some moderate multiple of £) and do not depend on the absolute 
field a(x). 

The SGS flux terms can be expressed in terms of increments by noting that gradient fields and central 
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moments are related to increments ai_ 

V7, = 0[5f{i)/i], 

n = o[5f{e)], 

Mf,g) = 0[5f{i) Sg{i)l 

VT,{f,g) = 0[Sf{i) 5gi£)/£], 

T,{f,g,h) = 0[5f{(.)5g{i)5h{(.)], 

where the symbol O stands for "same order-of-magnitude as" and 



(19) 
(20) 
(21) 
(22) 
(23) 



T'H/. ff. h) = {fgh)^ - fiTi{g, h) - geTe{f, h) - hiTi{f, g) - figghi (24) 

is the 3'"''-order generalized central moment of fields /(x), ^(x), /i(x) (see !25]). There are rigorous 
versions of relations p^ - (l^^ . See [H] and Appendix El below for details. 
Using relations (|19p and (PT|) . we can express baropycnal work as 



(25) 



In order to express 11^ in terms of increments, we write down the following identities which are 
straightforward to verify: 



u = u + r(p,u)/p 
djUi = djUi + p^^djT{p, Ui) - p~'^t{p, Ui)djp 
i,Uj) = T{ui,u.j) +^^^t{p,Ui,Uj) -^^'^t{p,Ui)t{p,Uj) 



(26) 
(27) 
(28) 



It then follows that we can express deformation work p2|) as 

Ilf (x) = - p [djUi + p~^d-jT{p, Ui) - p'^t{p, Ui)djp\ 

X [r(ui,Uj) + p"^r(p, Ui,Uj) - p"^T(p, Ui)r(p, Uj)] . 
From and relations ([T^ . (PT|) - (P^ . we can express deformation work as 

'6u{e) Sp{£) 6u{£) Sp^{e) Su{e) 



(29) 



p 



P 



Spie) 



p 



SpHe) 



e 



p p' 

Expressions ((25)) and pOI) are not heuristic estimates, but are based on rigorous versions of ([T^- 
whose details can be found in [19] and in Appendix |^ 



(30) 



^ Relation I I23I I is based on an unpublished exact expression due to G. L. Eyink (see 1201 ). We repeat the details in 
Appendix [A] 
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4 Locality of increments 



Since Hi and Ag can be expressed in terms of velocity, pressure, and density increments, it thus becomes 
sufficient to show that these increments themselves are scale-local. To establish this, we need the second 
requirement crucial for locality — that scaling properties of structure functions of velocity, pressure, and 
density increments are constrained by: 

< ct;; < 1, (31) 
< < 1, (32) 
< ct; (33) 

for some dimensionless constants Ap, Bp, and Cp. The root-mean-square of a field /(x) is denoted by 
frms = {f^y^^, where (...) is a space average, J dx{. . .). Here, an Lp-norm || • \\p — (| • |p)1/p is just the 
traditional structure function = (| • |^) raised to the 1/p-th power. The scaling conditions (PT|) - ([55)) 
are well-established empirically in incompressible turbulent flows over intermediate scales L r ^ 
They also have strong empirical support from many independent astronomical and numerical studies of 
compressible turbulent flows, which we discuss in detail in § 17.11 Note that condition (p3| on the scaling 
of density increments is only an upper bound. It only stipulates that the intensity of density fluctuations 
decays at smaller scales, which is a very mild requirement and is readily satisfied in incompressible or 
nearly-incompressible flows. 

Scaling conditions (PT|) - (|55)) reflect a structure of fields at intermediate scales. The constraints ap < 1 
and (7p < 1 indicate that velocity and pressure fields should be "rough enough" in a mean sense. They 
are not satisfied, for example, in laminar flows. These conditions are used to prove that contributions 
from the very large scales A to the flux across ^ ^ A are negligible or, in other words, that the flux is 
infrared local. 

On the other hand, constraints dp > and > indicate that the velocity and density fields 
are "smooth enough" in a mean sense. They are not satisfied, for example, if the fields are dominated 
by small-scale fluctuations with a non-decaying spectrum. These conditions are used to prove that 
contributions from the very small scales S to the flux across t ^ 5 are negligible or, in other words, that 
the flux is ultraviolet local. 

Notice that, unlike for the velocity and pressure fields, we do not stipulate that the density field be 
"rough enough." Contributions to the flux across scale ^ from the largest density scales L ^ £ need not 
be negligible, yet the flux can still be scale-local. The underlying physical reason is simple; an energy flux 
across scale £ at a point x will depend on the mass in a ball of radius i around x. Mass is proportional 
to average density ^^(x) in the ball, which is dominated by large scales: ^^(x) = ©[^^(x)] = 0[prms\- 



||(5u(r)||p - UrmsAp{r/LYp , 
\\5P[v)\\p ^ P„nsBp{r/L)< , 
||<5p(r)||p < Pr^sCpir/Lyi, 
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Indeed, for incompressible turbulence, the only density scale present is a k = Fourier mode, the largest 
possible, and the scale-local SGS flux is directly proportional to this density mode. 

Furthermore, we do not require that the pressure field be "smooth enough," even though we expect 
ap > 0. This is because pressure only appears as a large-scale pressure-gradient in (II3I) . For any filter 
kernel G(k) which is compact in Fourier space, contributions from wavenumbers Q 3> i^^ will be exactly 
zero and ultraviolet locality of VP^ is guaranteed without any scaling assumptions. 

Under conditions pip - ([55|) . proving scale locality of increments becomes simple and follows directly 
from 19J. For example, the contribution to any increment 6f{£) from scales A > £ is represented by 
Sf^{£). Here, /(x) can denote either velocity or pressure field. Since the low-pass filtered field /a(x) is 
smooth, its increment may be estimated by Taylor expansion, and (jl9p . and pip or p2p. as 



\\Sf^{£)\\p^\\£.{Vf 



A)\\P 



= o 



-(- 

A I L 



= O 



(34) 



and this is negligible for A 3> ^ as long as cr^ < 1. The notation 0{. . . ) denotes a big-O upper bound. 

On the other hand, the contribution to any increment Sf{£) from scales d < £ is represented by 
Sfg{£). Here, /(x) can denote either velocity or density field. Since fg = 0[f{x + S) — f{x)] from ([20]) . 
scaling conditions (|3T|) and (p3|) imply that 









= 





(35) 



and this is negligible for S £ as long as cr| > 0. For more details and for the careful proofs of these 
statements, see |19) . 

5 Locality of flux 

Based on results in §[3] and §[4l proving scale locality of the fiux terms Hf and is straightforward. 
To illustrate, consider fiux due to baropycnal work, Ai in (|13p . It is a quartic quantity which depends 
on two density modes, one pressure mode, and one velocity mode. This dependence can be made more 
explicit by writing 

Ae{p,P,p,u} = ^VPi-Ti{p,u), 

Pe 

where the first density argument A(p, ., ., .) corresponds to the factor 1/p. 

To prove ultraviolet locality of A^, we need to show that contributions from each of the four arguments 
(p, P, p, u) at scales S £ is negligible. It is obvious that will have vanishing contribution (decaying 
faster than any power, or exactly zero for filter kernels G'(k) compact in Fourier space) from scales S <ti £. 
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It follows that its reciprocal also has vanishinj^ contribution from very small scales 5 <^ I. It is 
also obvious that VP^ will have vanishing contribution from scales 6 <^ I. 

What remains to be shown is that the last two arguments u) have negligible ultraviolet con- 

tributions. If we replace them by p'^ and u^, respectively, we get from using Holder's inequality that 

||A,(p,F,p^,ui)||p < (const.)||^||oo||VP,||3pPp^WI|3p||'5iiiWll3p 



= O 



(36) 



which vanishes asS/£ for any Lp-norm if cr^p+cr^p > 0. We used relation (PT|) to get \\Ti{p'g, u'g) ||3p/2 < 
(const. )||J/3^(£)||3p||(5u^(£)||3p and result (I55t to arrive at the upper bound. A slight complication in our 
result is the additional overall factor {£/L)°', where a = cr^p + cr^p + cr^p — 1- This can grow with decreas- 
ing ^ if a < 0, causing our upper bound p6p to deteriorate at small scales. Such a growth, however, 
can be easily offset by taking 6 small enough: 6 < S^{e) = £ (^/i)(i"'^S>~'^3p-'^3p)/('^Sp+'^3p). We finally 
conclude that flux term is ultraviolet local under scaling conditions (PT|) . p3p for the velocity and 
density fields. 

Next we prove infrared locality of Ai (p, P, p, u) by replacing the pressure and velocity arguments 
with Pa and ua, respectively, for A ^ £. Using Holder's inequality, we get 

||Af(p,PA,P,UA)|lp < (const.)||^||oo||<5PA(^)/^||3p||M^)l|3p||'5uA(^)||3p 



= O 



L) V A 



(37) 



which vanishes as £/A for any Lp-norm if ct^ -|- a^p < 2. We used relations (fT9)) and pil) to get 
the inequality and result p4p to arrive at the upper bound. As in the ultraviolet case, there is a factor 
[t/L]'^ which can cause our upper bound (1571) to increase at small if a < 0. Again, such an increase 
can be easily compensated by taking A large enough: A > A,(£) = £ {^ij L)i'^ip+"3p+'^3p-^)/<~'^-'^3p-'^3p) ^ 
We finally conclude that fiux term Kg is infrared local under scaling conditions (|3T]) . ([32]) for the velocity 
and pressure fields. 

Similarly, we can derive rigorous upper bounds on non-local contributions for each of the 9 terms in 
p9| to prove scale locality of flux due to deformation work, H^. 



6 Locality of the cascade 

In proving locality of the SGS flux, H^ + A^, we did not need to make any assumptions about an equation 
of state for the fluid. We also did not analyze the internal energy budget ([3]). It is not obvious to us how 

^ For any positive smooth filter kernel, the field Pf{x) is real analytic and non-zero for all x. A well-known result 
from real analysis (see for e.g. 1330 states that the reciprocal is also real analytic and, therefore, its Fourier mode 

amplitudes decay faster than any power n of wavenumber fe~" as A: — >■ oo. 
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to best define a notion of scale for internal energy consistent with our scale-decomposition of the flow 
field in [1. Despite this shortcoming, we were still able to derive important results concerning locality. 

Circumventing the aforementioned shortcoming was possible due to two facts. First, we proved 
rigorously in [T] that viscous dissipation is negligible at large-scales, which implies that large-scale kinetic 
energy does not couple to internal energy via viscous dynamics. The only coupling that exists is via 
large-scale pressure dilatation. The second fact which aided us in proving locality of the SGS flux is that 
large-scale pressure dilatation, — P^V-u^, does not involve scales < t and, therefore, it cannot transfer 
kinetic energy directly across scales. Hence, pressure dilatation does not contribute to the SGS kinetic 
energy flux. 

In principle, one could conjure a scenario in which mean large-scale kinetic energy is converted to 
internal energy at scale ^ £i via PD{£) = — (P^V-u^) and is subsequently re-converted back, indirectly, 
into mean kinetic energy at a much smaller scale £2 ^ £i- In other words, as we continuously probe 
smaller scales £, the function PD{£) is positive at ^ £1, then decreases and becomes zero at ^ £2- The 
problem arises if such a process repeats itself, whereby PD{£) oscillates indefinitely with a non-decaying 
amplitude, as ^ — > 0. 

As unlikely as such a scenario may appear, we do not know of a rigorous argument to disprove it 
under the weak assumptions ()31|) -(|33 |) we have already made. It, therefore, precludes us from claiming 
at this point in the paper that our proof of a scale-local SGS fiux implies rigorously a scale-local cascade 
process in compressible turbulence. It is possible, however, to infer rigorously that the cascade is scale- 
local if we make one additional assumption which is, albeit reasonable, not as weak as scaling conditions 

dSlli-ia]). 

6.1 Pressure dilatation co-spectrum 

The assumption we need to prove a scale-local cascade concerns the pressure dilatation co-spectrum, 
defined as 

E^^ik)^ E -P(k)V:u(-k), (38) 

fc-0.5<|k|<fc+0.5 

which we require to decay at a fast enough rate, 

|P^^(fc)| < CUr^s Prras {kL)-^, /3 > 1 . (39) 

Here, C is a diniensionless constant and L is an integral scale. 

In the limit of infinite Reynolds number, assumption ([39]) implies that mean pressure dilatation, 
PD{£), asymptotes to a finite constant, 9 = — (PV-u), as £ — >■ 0. In other words, mean pressure 
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dilatation, PD{£), acting at scales > i converges and becomes independent off at small enough scales: 

lim PD{£) = lim V E^^{k) = 9, (40) 

for wavenumber K i^^ . In i) l7.2[ we shall give a physical argument on why we expect pO|) to hold. 

We remark that condition is sufficient but not necessary for the convergence of PD{£) in the 
limit of ^ 0. It is possible for E-^^{k), which is not sign-definite, to oscillate around as a function 
of k, such that the series J2k<K E^^{k) converges with /-C — > oo at a rate faster than what is implied 
by assumption (|39)). 



6.2 Conservative kinetic energy cascade 

Saturation of mean pressure dilatation, as expressed in (1401) . reveals that its role is to exchange large- 
scale mean kinetic and internal energy over a transitional "conversion" scale-range. At smaller scales 
beyond the conversion range, mean kinetic and internal energy budgets statistically decouple. In other 
words, taking i!^ ^> first, then £ ^ 0, the steady-state mean kinetic energy budget becomes, 

(n, + A,) = (£"'■'■) -0. (41) 

We stress that such a decoupling is statistical and does not imply that small scales evolve according to 
incompressible dynamics. However, while small-scale compression and rarefaction can still take place 
pointwise, they yield a vanishing contribution to the space- average. 

We denote the largest scale at which such statistical decoupling occurs by ic- It may be defined, for 
instance, as the scale at which PD{£c) — 0.95 9. Alternatively, it may be defined as 

_Ekfc"'^^^(k) 

= Ek£^^^(k) ■ ^^^^ 

Over the ensuing scale-range, £c > £ ^ £p,, net pressure dilatation does not play a role, and if, further- 
more, (e™^) in (|4T|) is localized to the largest scales as shown in ^j, then (Ilf -I- A^) wiU be a constant, 
independent of scale I. 

A constant SGS flux implies that mean kinetic energy cascades conservatively to smaller scales, 
despite not being an invariant of the governing dynamics. This is one of the major conclusions of our 
paper. In particular, the scenario discussed in the two paragraphs preceding ij lG.ll cannot ocure over 
tc > £ ^ £f_i, and kinetic energy can only reach dissipation scales via the SGS flux, Hf -\- A^, through 
a scale-local cascade process. We are therefore justified in calling scale-range £c > £ ^ £fj, the inertial 
range of compressible turbulence. 
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7 Discussion 



In proving locality of the SGS flux, 11^ + A^, we did not invoke assumptions of homogeneity or isotropy. 
We only assumed the weak scaling conditions pip - p3p on structure functions of velocity, pressure, and 
density. The results also apply to individual realizations of the flow, without the need for ensemble 
averaging. 

7.1 Validity of assumptions (l3B- (l33|) 

We have proved through an exact analysis of the fluid equations that scaling assumptions (I5T|) - (I551) are 
sufficient to guarantee scale locality of SGS flux. The ultimate source of these scaling properties is 
empirical evidence from experiments, astronomical observations, and numerical simulations 

For incompressible turbulence, which may be viewed as a limiting case of our analysis, the scaling of 
velocity and pressure structure functions has been well-established by a variety of independent 

studies such as those by 01 [5J [T^l HHl [23 HSl HI] • Assumption on density structure functions is 
trivially satisfied for a uniform density field. 

As for compressible turbulence, the available data is also in compelling support of our assumptions. 
Astronomical observations by [5J [B] of radio wave scintillation in the interstellar medium, characterized 
by highly supersonic turbulent flows, possibly up to Mach 20 (see for example [H]), show that 2nd-order 
density structure function scales with = 0-3 over 5 decades in scale. [50l[7] used velocity integrated 
spectral line maps of several molecular clouds and found evident power-law scaling for the density with 
exponent 0.3 < < 0.4. Much effort has also been expended to measure scaling of 2nd-order velocity 
structure functions. Using spectroscopic surveys of molecular clouds, which give line-of-sight velocities 
from emission lines, several independent studies by []Falgaroneetal92,MieschBally94,BruntHeyer02; [41] 
found scaling exponents = 0.4. More recently, [5S] measured scaling of structure functions up to 
6th-order and found cr^ = 0.54, cr| = 0.51, cr^ = 0.49, = 0.47, cr^ = 0.46, = 0.45. Analysis of 
solar wind data has also yielded < (7^ < 1, for 1 < p < 6 (see for example []Podestaetal07 and |46j). 

Alongside observational evidence, [33] US] HT] [131 [JS] carried out several independent numerical 
studies of forced compressible turbulence and calculated scaling of structure functions. Their simulations 
employed an isothermal equation of state and reached a range of high turbulent Mach numbers (based 
on Urrns of vclocity fluctuations), Mt — 2.5 — 10. They report power-law scaling exponents well within 
our required constraints, < and < (7^ < 1 for 1 < p < 6. One caveat of such simulations is 
that they do not resolve viscous dynamics explicitly but rely on numerical schemes to deal with shocks 
and dissipation. Based on such considerations, they may not be deemed direct numerical simulations 

^ Our scaling conditions Il31l l-| [33ll do not distinguish between compressive u"^ and solenoidal components of the 
velocity field. Whereas the contribution to structure functions in the incompressible limit will be predominantly from u", 
contributions from u'^ may be significant in general, as in the case of Burger's turbulence. 
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but rather implicit LES with an uncontrolled subgrid model and, therefore, are not as reliable. DNS 
of compressible turbulence cannot achieve simultaneously high Reynolds and Mach numbers due to 
resolution limitations. The largest DNS to date we are aware of is that by [33|- It was on a 1024'^ grid, 
had a Taylor Reynolds number Re\ — 300 and a turbulent Mach number Mt — 0.3 — still in the subsonic 
regime. 

We can also gain useful insights into the weakness of conditions (PT|) - (I55|) through exact mathematical 
considerations. From the definition of a structure function, = {\5f{£)\P), for any field /(x) and its 

scaling exponent, = hminff^o[ln'S'/(^)/ln(^/L)] (see [HI HTJ [H] ) , it is known that C/ is a concave 
function of p G [0, oo) (see for example [24l [20] for details). It follows that our scaling exponents 
ajj — C//p sre non- increasing functions of p (see |20j). For example, if we have cti < 1, then we are 
guaranteed Up < 1 for any p > I. Similarly, if Cg > for some g > 1, then Cp > for any p < q. Another 
known result states that if the pth-order moment of /(x) exists, (|/|'') < oo for p> I, then the pth-ordcr 
scaling exponent is non- negative, > (see [^ [20] ). 

To further put our scaling assumptions into perspective, a 2nd-order structure function 5*2 (^) ^ ^^'^^ 
is related to the spectrum E'"'{k) ^ fc^" with n = 2(72 + 1- Therefore, condition ctI > is equivalent 
to a constraint on the spectral exponent n > 1. This condition on n is the same as that required for 
a stationary velocity field in a bounded domain to have finite variance, (|up) < oo, and a power-law 
spectrum fc^" for /c G [fco, oo) in the limit Re oo (see [24] ) . 

7.2 Validity of assumption ( [59]) 

The existence of a scale-local cascade over an inertial range is the main conclusion of this paper. However, 
to reach our result, we made the important assumption (|39p which deserves more careful examination. 
Needless to say, the scaling of pressure dilatation co-spectrum is easily measurable from numerical 
simulations. The only reported measurement of this quantity we are aware of is by []Leeetal06 shown 
in their figure 6(b). The authors had the same purpose in mind; to check the scales at which pressure 
dilatation exchanges kinetic and internal energy. From their plot, they concluded that such an exchange 
takes place only at the largest scales. While their conclusion is in support of our postulate, their plot is 
on a log-linear scale which precludes the inference of such a conclusion. It is possible for the co-spectrum 
to scale with /3 < 1 in ([39]) while seeming to have most of its contribution from the largest scales. 
The point we want to emphasize here is that it is PD{£), the integral of the co-spectrum, which needs 
converge and have most of its contribution from the largest scales. We note that our condition does 
not require a power-law scaling — only that E^^{k) decays at a rate faster than ^ k^^ . 

It is not at aU trivial why one should expect PD{£) = — (P^V-u^) to converge at small scales. How 
can this be reconciled with the expectation that compression, as quantified by V-u, would get more 
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intense at smaller scales? Indeed, it has been observed numerically that (V'u)rms is an increasing 
function of Reynolds number (see for example []Leeetal91). The key point here is that our assumption 
(1391) concerns spatially averaged pressure dilatation. It is true that V'u(x), being a gradient, derives 
most of its contribution from the smallest scales in the flow. However, since PV-u is not sign-definite, 
major cancellations can occur when space-averaging. The situation is very similar to helicity co-spectrum 
in incompressible turbulence, H{k) = ^^. q g u(k)'Vxu(— k). Here, the pointwise vorticity, 

tL'(x) = Vxu, can also become unbounded in the limit of infinite Reynolds number. However, numerical 
evidence shows that (u-lo) remains finite with Reynolds number and the co-spectrum H{k) decays at a 
rate ~ fc^", with n « 5/3 > 1 (see for example [11] and []Kurienetal04). 

We can offer a physical argument on why PD{£) is expected to converge for £ ^ as a result of 
cancellations from space-averaging. The origin of such cancellations can be heuristically explained using 
decorrelation effects very similar to those studied in [21] and [2]. While the pressure in (PV-u) derives 
most of its contribution from the largest scales, V-u is dominated by the smallest scales in the fiow. 
Therefore, pressure varies slowly in space, primarily at scales ^ L, while V-u varies much more rapidly, 
primarily at scales ^ P, leading to a decorrelation between the two factors. More precisely, the 
pressure P^ in PD{£) may be approximated by P( — 0[Prms] — C>[Pl] such that 

(PiV-u,) « (PlV- ((n^)^ + (u,)^)) « (PlV-ul) + (Pl)(V.(u,)'^). 

The first term in the last expression follows from (u^)^ ~ u^, while the second term is due to an 
approximate statistical independence between Pl and V-(u^)'2^ ^ Su(£)/£ which varies primarily at 
much smaller scales i <ti L. If there is no transport beyond the domain boundaries or if the flow 
is either statistically homogeneous or isotropic, we get (V'(u£)'^) = 0. The heuristic argument finally 
yields that pressure dilatation, 

PP(£) = -(P,V-u,)«-(PiV-Ui), (43) 

becomes independent of £, for £ <^ L. Expression P5)) corroborates our claim that the primary role of 
mean pressure dilatation is conversion of large-scale kinetic energy into internal energy and does not 
participate in the cascade dynamics beyond a transitional "conversion" scale-range. 

7.3 A related study 

An insightful and clever numerical study by ;36i , which we mentioned in ^ \7.2l came to our attention 
at an advanced stage of writing this paper. The authors carried out DNS of compressible isotropic 
turbulence at low Mach number decaying under the influence of a randomly distributed temperature 
field. Some of the main conclusions are elegantly summarized by a schematic in their figure 8. They 
assert that mean pressure dilatation only acts at the largest scales beyond which mean kinetic energy 
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cascades conservatively down to the viscous scales where it is dissipated into heat. This conclusion is 
identical to the picture we arrived at in ii l6.2l and §[721 However, their study does not address the issue 
of scale locality which forms a main theme of our paper. 

On the other hand, their work goes beyond these statements and maintains that mean pressure 
dilatation couples internal energy to irrotational (and not solenoidal) modes of the velocity field. They 
also observe that mean pressure dilatation is not sign-definite in time but tends to transfer energy 
from internal to kinetic energy after time-averaging. Furthermore, they contend that the coupling of 
solenoidal and irrotational modes of the velocity field is weak and that each cascades separately to the 
viscous scales. They also investigate the alignments between vorticity and gradients of pressure, density, 
and temperature. All of these are essential issues which we do not tackle in our paper. 

The paper of [36) is a very valuable numerical investigation of the fundamental physics of compressible 
turbulence. Yet, inspection of dissipation spectra in their figures 6(c,d) seem to indicate that the 
turbulent flows they studied were not fully developed. We believe that similar studies at higher Reynolds 
and Mach numbers, and under different controlled conditions are still needed to establish the findings 
of [35] as empirical facts. 

7.4 Intermittency, locality, and universality 

As we mentioned above, scale locality is necessary to justify the notion of universality of intertial- 
range statistics. Kolmogorov's original 1941 theory of incompressible turbulence assumed statistical 
self-similarity of inertial-range scales. Today, there is a general consensus based on substantial em- 
pirical evidence that fiuctuations in turbulent fiows are not statistically self-similar but are subject to 
intermittency corrections. For p-th order structure functions, this is expressed as 

sM ^ {\su{i)n - (^^T ^ {^T^ ' ^^^^ 

where average energy fiux (e) (or dissipation) is empirically related to Urms and L through the zeroth- 
law of incompressible turbulence, (e) — u^^^^/L, and we have used Cp — p/3 + 5(p. Relation PH) shows 
that for non-zero "anomalous exponents" 6(^p, statistical averages at inertial-range scales £ L are a 
function of integral length, L. Colloquially, this means that inertial-range scales "remember" the number 
of "cascade steps", log2(i/^), for energy in going from scale L to scale £. 

However, intermittency does not contradict scale locality or the existence of universal scaling. It is 
known, for example, that the GOY shell model, in which scale interactions are local by construction, 
exhibits intermittency (see for example [9]). Despite remembering the number of cascade steps, there is 
no direct communication between inertial and integral length-scales due to scale locality. It is precisely 
because individual cascade step are scale-local and depend only upon inertial-range dynamics, that it 
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is possible to argue for universality of the scaling exponents Cp, regardless of (and consistently with) 
intcrmittcncy corrections. For a more detailed discussion of such issues, see [5D]. 

7.5 Shocks and locality 

An idea especially common in the astrophysical literature claims that in compressible turbulence a "finite 
portion" of energy at a given scale must be dissipated directly into heat via shocks rather than cascading 
in a local fashion (see for example [38]). Our analysis shows that large-scale kinetic energy can only 
reach dissipation scales through SGS flux, + , which we have have proved to be scale-local provided 
the weak scaling conditions (|3ip - ([55|) are satisfled. Therefore, in order for large-scale kinetic energy to 
dissipate into heat non-locally, it is necessary to violate ([5T |) . ([55)1 such that CTp < or (t|^ < for p < 7 
to break down ultraviolet locality. Having ap''' = implies that the mean intensity of velocity or density 
fluctuations does not decay at smaller scales. All empirical evidence discussed in subsection 17.11 seems 
to rule out such a possibility. 

The situation in compressible turbulence is similar to that of incompressible MHD turbulence where 
discontinuities in the magnetic field, i.e. current sheets, are pervasive. However, [3] proved rigorously, 
under scaling conditions analogous to pil) - ([55|) . that the cascade is local in scale and, furthermore, 
provided numerical support from high-resolution DNS. 

Another elucidating example is that of Burger's turbulence, in which viscous dissipation takes place 
only in shocks (which have zero volume when the Reynolds number is infinite) and vanishes everywhere 
else. However, scaling exponents of velocity increments are known to be cTp = 1/p for p > 1, which 
satisfy the condition < (7^ < 1 for any 1 < p < oo, and the proof of locality applies to Burger's 
energy cascade as well. In fact, the same conclusion was pointed out by [32] where he stated that the 
cascade in turbulence can be scale-local despite the presence of coherent discontinuous structures. He 
specifically discussed Burger's turbulence and showed that the energy is transfered to smaller scales in 
a local cascade process. As we mentioned in the introduction §[1] Kraichnan at the time had realized 
through his closure theory that scale locality depends on the exponent of the spectrum power-law and 
not on coherence properties. Furthermore, scale locality is perfectly consistent with the possibility that 
all dissipation takes place only in shocks and singular structures, over a set of zero volume in the limit 
oi fi ^ 0. It is well-known for this to be the case in Burger's fiow. The point here being that scale 
locality is an inertial-range property of SGS flux which transfers energy across scales and not of viscous 
dissipation which converts energy into heat. 
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8 Summary 



We proved that in compressible turbulence, the SGS flux responsible for direct transfer of kinetic energy 
across scales is dominated by scale-local interactions. This was achieved through rigorous upper bounds 
on the non-local contributions, and under weak assumptions (|3ip - (|33p . No assumptions of homogeneity 
or isotropy were invoked and the results hold for any equation of state of the fluid. 

We also showed that, in the limit of high Reynolds number, scale locality of the SGS flux implies 
a scale-local cascade of kinetic energy if condition on pressure dilatation co-spectrum is satisfied. 
In particular, condition p9p implies that beyond a transitional "conversion" scale-range, there exists an 
inertial scale-range over which mean kinetic and internal energy budgets statistically decouple and the 
mean SGS flux, + A^), becomes a constant, independent of scale £. Our result in §[S] demonstrates 
that kinetic energy cascades conservatively despite not being an invariant. 

We remark that the extent of the conversion range could be an increasing function of Mach number 
and/or the ratio of compressive-to-solenoidal kinetic energy. If so, it would have an immediate bearing 
on the interpretation of results such as those in figure 2 of [ID], where the Mach number is varied while 
maintaining a fixed Reynolds number. In such studies, an increasing Mach number may lead to the 
conversion range eroding away the finite inertial range present in simulations. Measurements of power- 
law exponents over the conversion range may not reflect an asymptotic scaling which would otherwise 
appear at sufficiently high Reynolds numbers. The problem is that of an ordering of limits; the physically 
interesting order being one in which we take Re oo first, followed by Mt — > oo. 

In summary, we conclude that there exists an inertial range in high Reynolds number compressible 
turbulence over which kinetic energy reaches dissipation scales through a conservative and scale-local 
cascade process. This precludes the possibility for transfer of kinetic energy from the large-scales directly 
to dissipation scales, such as into shocks, at arbitrarily high Reynolds numbers as is commonly believed. 
Our work makes several assumptions and predictions which can be tested numerically. Our locality 
results concerning the SGS flux can be verified in a manner very similar to what was done in [2] and 
[3]. We also invite numerical tests of assumption p9p on the scaling of pressure dilatation co-spectrum. 
Verifying ((39| or (|40)) under a variety of controlled conditions would substantiate the idea of statistical 
decoupling between mean kinetic and internal energy budgets. This would have potentially significant 
implications on devising reduced models of compressible turbulence as well as providing physical insight 
into this rich problem. 
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A Details of relations (IT9|)-(l23|) 



We repeat from [TH] and [20] the details of expressing gradients and central moments in terms of incre- 
ments. 

Relation ([T9| is a result of integration by parts and rewriting a large-scale gradient as: 

V7,(x) = -jjdv{VGUv)f{^ + v) 

= -iy'dr(VG),(r)[/(x + r)-/(x)], (45) 

where {'VG)i{r) — £~''(VG)(r/i?) in d-dimensions and we have used condition / c!rVG(r) = to arrive 
at the last equality. 

High-pass filtered fields in relation ((20)) can be expressed in terms of increments as: 



/;(x) = -y drG,(r)[/(x + r)-/(x)] 

= -(<5/(x;r)},, (46) 

where ((5/(x;r))f = J drG£{r)df{x;r) is an average over separations r in a ball of radius of order i 
centered at x. 

The 2"'*-order central moment relation (|2ip is due to [] Constantinetal94 and is straightforward to 
verify: 

T, (/ (x) , g (x)) = (J/(x; r)59(x; r)), - (,5/(x; r)),{Sg{^- r)),. (47) 
Relation ([2^ follows from (|T7)) through an integration by parts: 

VT,(/(x),.g(x))--i| I dr(VG),(r)<5/(x;r)59(x;r) 

- J dr(VG),(r)<5/(x;r) J dr'G,(r')<5.9(x; r') 

-J drG,(r)5/(x;r) J dr'(VG),(r')<5.9(x; r')|. (48) 
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The 3'''*-order central moment in relation (1^^ is due to f5(T and is straightforward to verify: 



((5/(x; r)Sg{x; r)6h{x; r))e 



(<5/(x;r)),(5g(x;r)<5/i(x;r)), 



(<55(x;r)),(5/(x;r)5/i(x;r)}f 



(<5/i(x;r)),((5/(x;r)(5g(x;r)), 



+ 



2(J/(x;r)),(<5g(x;r)),(J/i(x;r)),. 



(49) 



B Locality proofs 

The proofs for locality follow closely those of HKi]. We will repeat briefly the ones which were explicitly 
discussed in the latter work and also give some details on locality of new terms not present in the [19] 
treatment. 

The filter kernel G{r) used in the proofs is smooth and decays faster than any power r^P as r — 7> oo. 
ft is also possible (but not required) to have both it and its Fourier transform G'(k) be positive and, 
furthermore, to have the latter compactly supported inside a ball of radius 1 about the origin in Fourier 
space. 

B.l Ultraviolet locality 

Ultraviolet locality means that contributions to the energy flux across £ from scales S £ decay at least 
as fast as 6°", for some a > 0. We will now show that each of the factors in SGS flux terms (fT3|) and 
([29|) is ultraviolet local. 

Non-local ultraviolet contributions to a large-scale gradient of field /(x) can be shown to be bounded 
using p5|) - p6)) as was proved by [19] : 



in Fourier space. 

Non-local ultraviolet contributions to a second-order central moment of fields /(x) and g{x) can be 
bounded using (H5)) . (|T7)) as was proved by [12]. For l/p = 1/r + 1/s, we have 




(50) 



Notation 0{. . . ) denotes a big-O upper bound. In fact, || Vu'^||p — for a filter G'(k) compactly supported 




(51) 
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Similarly, using (P^ . (|48p we can show that non-local ultraviolet contributions to the gradient of a 
second-order central moment are bounded. For l/p = 1/r + 1/s, we have 

||Vr,(/^,g^)||, < 1( / dr|(VG),(r)| ||/Kx)||,||.gKx)|U 



dr|(VG),(r)| \\n{^)\\rj dr\G,{r)\ \\g's{^)\ 
+ J dr|G,(r)| \\n{^)\\rj rfr|(VG),(r)| ||5^(x)| 

O -^]- (52) 



The non-local ultraviolet contributions to a third-order central moment of fields /(x), ^(x), and h(x) 
can be bounded using (|46p .(H ^ . For l/p = 1/r + 1/s + 1/t, we have 



\Mf5,9's:h's)\\p < 8 j dv\G,{v)\ ||/Kx)||,||g^(x)|U||/i^-(x)||t 

+ 8 j dv\G,{v)\ m^)\\rj dr\Gdr)\ \\g'si^)UK{^)\\t 

+ 8 J dr|G,(r)| Ui^)\\r J dr\G,{r)\ ||/Kx)||,||/j^(x)||, 

+ 8 / dr\Gdr)\ \\K{^)\\r f dr|G,(r)| \\f^{^)\\Mi^)\\t 



+ 16 / dr|G,(r)| ||/Kx)||. / dr|G,(r)| / dr|G,(r)| ||/z^(x)||, 



= O ((5^'+-='+-'") . (53) 
B.2 Infrared locality 

Infrared locality means that contributions to the energy flux across £ from scales A ^ £ decay at least 
as fast as A^", for some a > 0. We will now show that each of the factors in SGS flux terms ([TS]) and 
([29| . except for the density, is infrared local. As we have discussed above, it is physically expected that 
the kinetic energy cascade will depend on density fluctuations at the largest scales. 

Non-local infrared contributions to a large-scale gradient of field /(x) can be bounded using P5|) as 
was shown by |19j : 

II(V7a)J1p< / dr'|G,(r')| ^ J dr|(VG)A(r)| ||5/(x; r)||p ^ O (a-.^-^) . (54) 
Non-local infrared contributions to a second-order central moment of fields /(x) and (7(x) can be bounded 
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using (|47l) as was shown by [T^. For 1/p = 1/r + 1/s, we have 

I|tK7a,5a)IIp < / dr|G,(r)| ||(57^(x; r) ||,||<5g^(x; r)| 



+ j dv\G,{v)\ \\5f^{^-v)\\r j dv\Gt{v)\ \\5g^{^-v)\l 

where the last step follows from ||(55^(x; r)||s = 0{r/S."'i~^) (see [H]). 

Using (|3S]), we can also show that non-local infrared contributions to the gradient of a second- 
order central moment are bounded. Such a term only appears in the form V-T£(p, u) in (j29l) . For 

1/p = 1/r + 1/s, we have 

||V-t,(p,ua)||p < dr|(VG),(r)|. ||<5p(x;r)||,PuA(x;r)|U 

dr|(VG),(r)| ||<5p(x;r)||.- J dv\Ge{r)\ ||<5ua(x; r)||, 

dr\Ge{r)\ \\6p{^;r)\\r J dr|(VG),(r)H|(5uA(x; r)|u| 

O (a<-i K) . (56) 

We can also show infrared locality of a third-order central moment. Such a term appears in deforma- 
tion work (|29p in the form Ti(p, u,u). Non-local infrared velocity contributions can be bounded using 
(|49l) . For 1/p = 1/r + 1/s + 1/t, we have 

||r,(AUA,UA)I|p < J dr|G,(r)| ||V(x; r)I|,||5uA(x; r)IUPuA(x; r)||t 

dr|G,(r)| ||<5p(x;r)||,.y dr|G,(r)| ||<5ua(x; r)|U||<5uA(x; r) ||t 
dr|G,(r)| ||<5uA(x;r)||,y dr|G,(r)| H^x; r)|U||<5uA(x; r)||i 
dr|G,(r)| ||5uA(x;r)||,y dr|G,(r)| H^x; r)||.||<5uA(x; r)||i 
+ 2 j dv\Gt{v)\ \\6p{^;r)\\rJ dr|G,(r)| ||,5ua(x; r) |j, ^ rfr|G,(r)| ||(5ua(x; r)||* 
= o(a<+<-2 . (57) 
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